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Abstract 

In this Letter, the renormalization-group equations for the (flavor-conserving) CP-violating interaction 
are derived up to the dimension six, including all the four-quark operators, at one-loop level. We apply 
them to the models with the neutral scalar boson or the color-octet scalar boson which have CP-violating 
Yukawa interactions with quarks, and discuss the neutron electric dipole moment in these models. 
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I. INTRODUCTION 



The electric dipole moment (EDM) for neutrons is sensitive to CP violation in physics beyond 
the standard model (SM) around TeV scale. This is because, while the CP phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix is O(l), the CKM contribution to the neutron EDM is too much 
suppressed |Q] to be observed in near future. (The recent evaluation of the CKM contribution to 
the neutron EDM is given in Refs. The naturalness problem in the Higgs-boson mass term in 
the SM might require new physics at TeV scale, and many extensions of the SM generically have 
CP-violating interactions. The supersymmetric standard model, which is the leading candidate for 
the TeV-scale physics, is severely constrained from the EDM measurements Q . 

The (flavor-conserving) CP-violating effective operators at parton level up to the dimension six 
are the QCD theta term, the EDMs and the chromoelectric dipole moments (CEDMs) of quarks, the 
Weinberg's three-gluon operator Q and the four-quark operators. In the evaluation of the neutron 
EDM, the CP-violating four-quark operators tend to be ignored since the four-light-quark operators 
suffer from chiral suppression in many models. However, the four-quark operators including heavier 
ones, such as bottom/top quarks, may give sizable contributions to the neutron EDM. The EDMs, 
CEDMs, and the three-gluon operator are radiatively generated from the four-quark operators by 
integrating out heavy quarks. 

In the multi-Higgs models, the Barr-Zee diagrams are known to give the sizable contribution 
to the neutron EDM |J. In the Barr-Zee diagrams the heavy-quark loops are connected to light- 
quark external lines by the neutral scalar boson exchange so that the CEDMs for light quarks are 
generated at two-loop level at 0(a s ). However, it is not clear which renormalization scale should 
be chosen for a s . In addition, the contributions from the Barr-Zee diagrams at two-loop level to 
the quark EDMs vanish at 0(a s ). However, it is still unclear that the higher-order corrections to 
the quark EDMs are negligible in the neutron EDM evaluation. 

In this Letter, in order to answer those questions, we derive the renormalization-group equations 
(RGEs) for the Wilson coefficients for the CP-violating effective operators up to the dimension six 
at one-loop level, including operator mixing. The RGEs for the EDMs and CEDMs for quarks 
and the three-gluon operator have been derived in Ref. The next-leading order corrections 

to them are also partially included ||. We include the four-quark operators in the calculation at 
the leading order. Using the derived RGEs, we evaluate the EDMs and CEDMs for light quarks 
and the three-gluon operators induced by the neutral scalar boson exchange including the QCD 
correction. We also discuss the four-quark operators induced by the color-octet scalar boson. 
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This Letter is organized as follows. In the next section, we review the neutron EDM evaluation 
from the parton-level effective Lagrangian at the hadron scale. In Section 3, we derive RGEs for 
the Wilson coefficients for the CP-violating effective operators up to the dimension six at one-loop 
level. In Section 4, we show the effect of the running a s on the evaluation of the Wilson coefficients, 
assuming the neutral scalar boson exchange induces the CP-violating effective operators. In Section 
5, another example of the application of the RGEs is shown, assuming the effective operators 
induced by a color-octet scalar boson. Section 6 is devoted to conclusion. 



II. NEUTRON EDMS 



First, we review about evaluations of the neutron EDM from the low-energy effective Lagrangian 
at parton level. The CP-violating interaction at parton level around the hadron scale (fi H = 1 GeV) 
is given by 



i-CPV = ^TT^uv^ 



a s 

2 Yl d Q^( F " cr )75<? - I d q qg s {G ■ a)j 5 q 

q=u,d,s q=u,d,s 



+ IwfABoG^G^G?. (1) 

Here, F^ v and G A V {A = 1-8) are the electromagnetic and gluon field strength tensors, g s is the 
strong coupling constant (a s = g^/iir), F-a = F^a^ , G-a = G A u a pu T A , and G A V = \e^ upa G Apa 
with = §[7^, 7"] and e 0123 = +1. The matrix T A denotes the generators in the SU(3)c 
algebra, and f ABC is the structure constant. The first, second, third and forth terms in Eq. (|l|) 
are called the QCD 9 term, the EDM and the CEDM for quarks, and the three-gluon operator, 
respectively. In this Letter, the covariant derivative is defined as = — ieQ q A^ — ig s G A T A , 
in which and G A are gauge fields for U(1)em and SU(3)c, respectively with Q q , the QED 
charge ({Q u ,Qd,Qs) = (2/3,-1/3,-1/3)). In Eq. (||), we ignore the CP-violating four-quark 
operators, since their coefficients are often proportional to the light-quark masses in typical models, 
as mentioned in the Introduction. 

The neutron EDM is evaluated from the low-energy interaction at parton level with the naive 
dimensional analysis, the chiral perturbation theory, and the QCD sum rules, though they are 
considered to have large uncertainties. The evaluation in term of the QCD sum rules is more 
systematic than the others, at least for the contributions from the QCD theta term, and the quark 



EDMs and CEDMs to the neutron EDM [10]. The recent evaluation of the neutron EDM with the 
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QCD sum rules Q is 

d n ~ 2.9 x 10~ 17 (9 [e cm] + 0.32d d - 0.08d u + e(+0.12d d - 0.12(4 - 0.0064) • (2) 

In the evaluation, the recent QCD lattice result is used for the low-energy constant A n , which is 
defined by {0\rj n (x)\N(p, s)) = X n u n (p,s) with r/ n (x) the neutron-interpolating field. If a value of 
\ n evaluated with the QCD sum rules is used, the neutron EDM is enhanced by about five times 
compared with Eq. (||). 

The contribution from the three-gluon operator might be comparable to the quark EDMs and 
CEDMs. The quark EDMs and CEDMs are proportional to the quark masses, while the three-gluon 
operator does not need to suffer from chirality suppression. However, the size of the contribution 



from the three-gluon operator depends on the methods of the evaluation. In Ref. [12| the authors 
compare the several evaluations and propose 

d n (w) ~(10 - 30)MeV x ew . (3) 



III. OPERATOR BASES AND ANOMALOUS DIMENSION MATRIX 

We would like to introduce heavy quarks in the low-energy effective theory and evaluate their 
contributions to the neutron EDM. In this section, we show the one-loop RGEs for the Wilson 
coefficients for the CP-violating effective operators up to the dimension six, including heavy quarks. 

First, we define the operator bases for the RGE analysis. The flavor-conserving effective oper- 
ators for the CP violation in QCD are given up to the dimension six as 

£cpv= E £c?( M )0?(m) + C' 3 (m)<W 

i=l, 2,4,5 q 

+ EE cf 9 M% q M + \ E E cf 9 w&f 9 w ■ w 

i=l,2 q'jiq i=3,iq'^q 

where the sum of q runs not only light quarks but also heavy ones, and we ignore the QCD theta 

term since it is irrelevant to our discussion here . The effective operators are defined as 

i 

0\ = - -m q qeQ q (F ■ a)-y 5 q , 

1 

°l = - 2 m ^9s(G ■ a)j 5 q , 

3 = - \g s f ABC e^G A ^G B v X G c pa , (5) 

1 The QCD theta term does not contribute to the RGEs for other CP-violating terms. Furthermore, there may 
be contribution to the QCD theta term from other CP violation terms, while the QCD theta term vanishes 
dynamically if the Peccei-Quinn symmetry is invoked. 
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and 



oi 


= q a q a qp ijsqp , 




— qa° qaH/3 lu (j,v 75y/3 


dp 


= q' a ( i'aW i 'y5q/3, 


op 


= q' a q'i3qp i i5q a , 


dp 


= q' a ^ v q'om ia ^qp 


op 


= q' a a^q'pqpio^q a 



(6) 

Here, m q are masses for quark q. In Eq. (@) we explicitly show the color indices, a and f3. A 
factor of 1/2 appears in front of the fourth term of Eq. (|j), since the term is symmetric under the 
exchange of q' and q. The Wilson coefficients in Eq. (^) are related to the parameters in Eq. ([!]) as 

d q = m q eQ q Cf(fi H ) , 



d q =m q C%(n H ) , 



w 



-g s C 3 (n H ) . 



The RGEs for the Wilson coefficients of these operators are given as follows, 

where the Wilson coefficients are written in a column vector as 

p _ (f^i r<i n ni ni ntfi ni'i r<<m' nil! r<<i'<i r" ' 

— [yi, 2 , ^3, U 4 , Ug, lyj , Lv 2 , , u 2 , >^3 , 

The anomalous dimension matrix is calculated at one-loop level as 





V) 







4tt7/. 



where 
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(13) 



and 



7s/ 



+4 +4 
-32iV - 16 -16 










1U m, Qq 















(14) 



(15) 



where C F {= (N 2 — 1)/(2N)) is the Casimir constant of the fundamental representation, N(= 3) 
is the number of the color, rif is the number of light flavor quarks, and /3o(= 11/3 x N — 2/3 x nj) 
is the leading-order beta function of strong coupling constant. 

The anomalous dimensions for the dimension- five operators are calculated in Ref. ||, and that 
for the three-gluon operator is calculated by Ref. jj]]. The mixings among the dimension- five 
operators and the three-gluon operator are found in Ref. ||. We newly calculate other terms in 
the anomalous dimensions matrix, which are related with the four-quark operators. We note that 
the operator mixings between the EDM or CEDM operators and the four-quark operators are 
generated at 0{oP s ). 

IV. NEUTRAL SCALAR BOSON EXCHANGE 



In multi-Higgs models, a color-singlet neutral scalar boson (j) may have the CP-violating Yukawa 
coupling with quarks. If the Yukawa interaction violates the CP invariance, the CP-violating four- 
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quark operators are induced at tree level, after integrating the neutral scalar boson out, as 

2 

C\ = V2Gf — |/|/p i 
mi ° 



m cj) 

Cf = V2G F ^^flr p , (16) 



where we assume that (p is heavier than heavy quarks (m^ 3> m q ,m q '). Here, fg and f p are the 
CP-even and odd Yukawa coupling constants, respectively, defined as 

£ = 2 l l i G l F l2 m q q^(f q s + i f q p -y 5 )q a </>, (17) 

where (ft is a (CP-even) real scalar field, and Gf is the Fermi constant. We parametrize the Yukawa 
coupling constants as they are proportional to the quark masses. In typical new physics models, 
these Yukawa coupling constants are taken to be proportional to masses of quarks in order to avoid 
the stringent experimental constraints from the flavor physics data. For the SM Higgs boson, the 
Yukawa coupling constants are of /| = 1 and f p = 0. For the multi-Higgs models, e.g., two-Higgs- 
doublet models, the coefficients fg/ P may be much larger than unity because of the enhancement 
factor originated from the ratio of vacuum expectation values. 

It is known that, in these models, the EDMs and CEDMs for light quarks are generated by the 
Barr-Zee diagrams at two-loop level, and the three-gluon operator is also induced by the heavy- 
quark loops at two-loop level. Let us derive those contributions using the RGEs for the Wilson 
coefficients. 

In the leading-logarithmic approximation, since C\ is non-zero, the coefficients for the EDM 
and the CEDM operators are generated from the one-loop RGEs as 

C\ = CI = ~Cl In ^. (18) 

These contributions are compared with the explicit calculation of one-loop diagrams with the scalar 
boson exchange as 



(19) 



where a limit of m q <C m^ is taken. The second term in the parentheses should be considered as 
the short-distance contribution in which the loop momentum is around m^. 

The one-loop contributions from C| to the EDMs and CEDMs for light quarks (q = u, d, s) 
are negligible in the neutron EDM since they are suppressed by powers of the light-quark masses. 



However, when the CEDMs for heavy quarks are generated, the three-gluon operator is induced 
by integration of heavy quarks as follows [Of . 



c 3 K) = £^cfK). (20) 



Thus, from Eqs. ([L9L 2C), we get 



Ci—^V?! ft (21) 
Hit 5 \ m q 4 J 

The result is consistent with the explicit calculation of the two-loop diagrams for the three-gluon 
operator [14]. 

When C\ q and/or C qq are non-zero, the contribution to the CEDMs for light quarks is derived 
at the two- loop level, using the RGEs for the the Wilson coefficients, Eq. (|l0|) , as 

C q = ^^(ln^Y 



87r 3 m q V rriqi ) 



C qq + C qq . (22) 
This is because Cf q and C qq are mixed with C\ q and C qq which induce the CEDMs. When the 



Yukawa coupling constants in Eq. (17) are proportional to the quark masses, the induced CEDMs 
for light quarks are not suppressed by their masses, compared with the one-loop contribution. 
The result in Eq. (|22|) is consistent with the explicit calculation of the Barr-Zee diagrams M in a 



limit of S> m q ,m q i. It is known that the contribution to the neutron EDM from the Barr-Zee 
diagram at two-loop can be competitive with those from the one-loop diagrams || . Therefore, 
the four-quark operator for heavy quarks can be important in some cases. 

On the other hand, the EDMs for light quarks are not generated by two-loop level diagrams 
of 0(a s ), even if C\ q and Cf are non-zero. The EDMs for light quarks have contributions from 
C| 9 (C^ q ) and C\ q (Cf 1 ) in Eq. (|To|). Their contributions are exactly canceled with each others 
so that the EDMs vanish at the order. This is also consistent with the explicit calculation of 
the Barr-Zee diagrams. The EDMs generated at two-loop level are suppressed by a. However, 
the running effect of the strong coupling constant may prevent the cancellation so that the EDM 
would be enhanced. 

Now, let us consider the effect of the running strong coupling constant, a s (fx). Here we compare 
values of the EDM and CEDM operators for down quark and the three-gluon operators including 
and not including the renormalization-group evolution of the strong coupling constant. We assume 
that the Yukawa coupling constants for down and bottom quarks with <f> are non-zero in Eq. ([171) 
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and then 

Ci(^) = C 2 6 (m^) = + JLcf(m^) . (23) 

The last assumption comes from Eq. (|l9|). 

In Fig. 1 the CEDM for down quark, dd, (a) and the coefficient of the three-gluon operators, 
w, (b) at the hadron scale (fi = fi H = 1 GeV) are shown as functions of with fs = fp = ^ an d 
fg = fp = l. Here, we ignore the contributions from top quark, and other short-distance effects. 
If the scalar mass is larger than the top quark mass (m^ > m t ), the RGEs are solved using (3q 
with n f = 6, or if not, with nj = 5. When bottom quark is integrated out, the Wilson coefficient 
of Weinberg operator emerges. Then the RGEs are solved using (3q with nj = 4 to the scale 
/i = m c , and with n/ = 3 to the scale fi = 1 GeV. We use md(nH) = 9 MeV, m c (m c ) = 1.27 GeV, 
mb(mb) = 4.25 GeV, m t {m t ) = 172.9 GeV, and a s (m z ) = 0.12. For the coefficient w, we multiply 
10 MeV in the figure, which is a factor in Eq. (^), so that one may estimate the contribution to 
the neutron EDM. It is from Eqs. found that the three-gluon operator might be comparable 
to the CEDM when fj /p ~ fL p . 

In Fig. 2 the ratios of the CEDM for down quark (a) and the three-gluon operator (b) at 
/i = nib between including the running effect of a s and not including it (using the constant coupling 
a s = a s (m b )), are shown as functions of m It is found that the running coupling a s (u.) changes 
the CEDM by about 20% while the three-gluon operator is changed by at most 10 %. These results 
come from inclusions of the four-quark operators to the RGEs for the Wilson coefficients. 

In Fig. 3 the ratio of the EDM and CEDM for down quark is presented as a function of m^. 
The non-zero value of the EDM is generated as we mentioned. It is also found that the ratio is 
roughly proportional to logm., and the absolute value is about 0.14. In the evaluation of the 
neutron EDM with the QCD sum rules, the size of the contribution from the down-quark EDM is 
30-40 % of that from the CEDM. 

V. COLOR-OCTET SCALAR BOSON EXCHANGE 

In the previous section it is shown that the neutral scalar boson does not generate sizable 
EDMs for light quarks via two-loop diagrams at 0{a s ). This comes from cancellation between 
contributions via C| q (C^ 9 ) and C\ 9 (C% q ). This situation is different when the scalar boson 
has color. Now let us assume that a color-octet scalar boson £(= T, A T A ) is introduced and it has 
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FIG. 2: (a) : Ratio of the CEDM for down quark, dd, at yu = m& between including and not including 
running of the strong coupling constant, as a function of m^. (b) The same ratio for coefficient of three- 
gluon operator, w. 



CP-violating Yukawa interactions with quarks, 



(24) 



The octet scalar fields may appear in new physics beyond the SM such as a radiative seesaw model 
for the neutrino masses [15], and a grand unified model pi 
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Integration of the £ leads to the following Wilson coefficients, 

of. = _V2 Gf ^ ^/f/j- , cf « = ^ Gf ^I/j'/^ . 

cf = -V2G F ^±f s 4 , Cf = Vte F ^\tif$ • (25) 

Now C^ 9 and/or C| 9 are generated. At the two-loop level, the EDMs for light quarks have 
contributions via C| 9 (C 99 ) and C\ 9 (Cf 1 ) with different weights. This is different from the case 
of the neutral scalar boson exchange as discussed in the previous section. The EDMs for light 
quarks are generated at two- loop level of 0(a s ) as 

C 9 = -^^^C F (ln^) 2 \c 9 ' 9 + cf] , (26) 



47T 3 Tflq Qg x „„ ? . 

which is compared with the CEDMs for light quarks as 



C| = ^^fln msN 
87r d m 



~ C! ' 9 + i ^N~ Cf)CI ' 9 " ^ + i JN~ Cf)CQ A ■ (27) 



In Fig. 4 the EDM and CEDM for down quark at the hadron scale are shown as functions of 
m s with fg = fp = l and /| = f P = l 2 . Here we use the running coupling for a s . It is found 
that the EDM contribution is larger than the CEDM ones in the neutron EDM when we adopt the 
QCD sum rule result on the neutron EDM evaluation. (See Eq. (Q).) 

For completeness, we also show the magnitude of the three-gluon operator in Fig. 5. Here, we 
ignore the short-distance contributions to the CEDMs for heavy quarks and three-gluon operator 



If the electrically charged color-octet scalar boson exists, we may have large contributions from the top quark Jn|. 
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whose loop momenta are around mj and take C|(m E ) = Cs(m s ) = for simplicity 3 . Again, it is 
found that the three-gluon operator might give a comparable effect to other contributions. 
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FIG. 5: Three-gluon operator w at hadron scale as a function of m E . 



VI. CONCLUSION 

In this Letter, we have derived the renormalization-group equations for the CP-violating inter- 
action including the quark EDMs and CEDMs and the Weinberg's three-gluon operator as well as 
all the flavor-conserving four fermion operators. The operator mixings between the (C)EDM oper- 
ators and the four-quark operators are arisen at the order of a^, which can give large contributions 

3 The short-distance contribution to the three-gluon operator could come from a diagram in which three gluons are 
emitted from each quark and scalar boson line, and the evaluation is beyond the scope of this work. 
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to the EDMs via the renormalization-group evolution. 

Assuming the CP-violating Yukawa interactions for the neutral scalar bosons, it is known that 
the CEDMs for light quarks are generated from the diagrams with heavy-quark loops, called as 
the Barr-Zee diagrams. We show that when the neutral scalar boson is much heavier than heavy 
quarks, the Barr-Zee diagrams are systematically evaluated with the RGEs of the CP-violating 
interaction. We also show that the running effect of the strong coupling constant gives corrections 
to the contribution with more than 20 % compared with assuming the constant coupling. The 



uncertainties in the calculation of the neutron EDM have been estimated in the literature [17]. 
It gives about 50 % error for the QCD sum rule, while 40 % error for the low-energy constant 
evaluated from the lattice QCD calculation. Therefore, hadronic uncertainties would overcome 
the QCD corrections from the renormalization-group evolution at this moment. We hope that the 



lattice QCD simulation will improve and reduce uncertainties significantly [18]. 

The Barr-Zee diagrams at two-loop level do not contribute to the EDMs for light quarks at 
0(a s ). We show using the RGEs of the CP-violating interaction of a color-singlet scalar boson with 
quarks that ratio of the quark EDM over the CEDM is about 0.14, and it is roughly proportional to 
\ogm,. Thus, the contribution is not negligible to the neutron EDM at all. When the color-octet 
scalar boson has CP-violating Yukawa interaction with quarks, the quark EDMs are generated at 
two-loop level, and they are comparable to the quark CEDMs. 
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